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Abstract 

I compute the complete two-loop effective potential for the minimal supersymmetric stan- 
dard model in the Landau gauge. This enables an accurate determination of the mini- 
mization conditions for the vacuum expectation values of the Higgs fields. Checks on the 
result follow from supersymmetric limits and from renormalization-scale invariance. The 
renormalization group equations for the field-independent vacuum energy and the vacuum 
expectation values are also presented. I provide numerical examples showing the improved 
accuracy and scale dependence obtained with the full two-loop effective potential. 
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1 Introduction 

The mechanism of electroweak symmetry breaking will be the principal focus of experimental 
investigations at the high-energy frontier for the next decade. Supersymmetry provides a highly 
predictive mechanism for addressing the hierarchy problem associated with the electroweak sym- 
metry breaking scale. If supersymmetry is correct, then interpretations of future experimental 
data will rely on precise theoretical calculations in candidate models of supersymmetry breaking. 
The effective potential |]]-|| approach allows the computation of the vacuum expectation values 
(VEVs) of Higgs fields in terms of the underlying Lagrangian parameters of a given theory. 

The scalar potential of the Minimal Supersymmetric Standard Model (MSSM; for reviews, 
see is notoriously sensitive to radiative corrections. At tree level, the Higgs field quartic 

couplings are proportional to a sum of squares of electroweak gauge couplings and are therefore 
known not to be very large. Furthermore, they actually vanish along a D-flat direction in field 
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space. These facts ensures the existence of at least one light Higgs scalar boson, corresponding 
to a shallow direction in the effective potential for the Higgs vacuum expectation values. The 
same facts also imply that the minimization conditions for the scalar potential depend very 
significantly on radiative corrections. 

Previous results for the effective potential in the MSSM have included the full one-loop con- 
tributions and partial two- loop corrections [§]-|[T| including the effects of the QCD coupling 



and the top and bottom Yukawa couplings. Including these contributions mitigates the scale- 
dependence of the tree- level effective potential. However, I find that there is still a significant 
scale dependence and overall error compared to the uncertainties in theoretical quantities that 
may eventually be obtained at future experiments, especially at a linear e + e~ collider. In this 
paper, I will present the result for the full two-loop effective potential of the MSSM, in the 



Landau gauge and in the DR [12| regularization and renormalization scheme. This is an appli- 



cation of the results given in ref. |L3| for a general field theory. The DR scheme is the variant 



of the DR scheme [[T4| in which the effects of unphysical epsilon-scalars masses are removed by 
parameter redefinitions [|T|, |13[| . 

The two-loop effective potential for a general renormalizable theory can be written as 



V, 



cff 



167T 2 (167T 2 ) 2 

Here is the tree-level contribution. In the DR' scheme, the one-loop contribution is 

V^=J2(-l) 2Sn (2s n + l)h(ml), 



;i-2) 



where s n = 0,1/2,1 for real scalars, two-component fermions, and vector fields, respectively, 
with field-dependent tree- level squared masses m 2 . The one-loop function is 

„2 



h(x) = ^ [Hx/Q 2 ) - 3/2 
where Q is the renormalization scale. The two-loop contribution always has the form 

, npq 1 2 



;i.3) 



V (2) = J2\ nnpp F np (m 2 n ,m 2 p ) + £ \X^rFn Pq [m- n ,m-,m-), 



;i-4) 



2 2 2\ 
n,p n,p,q 

where A npi?r and \ npq are tree-level field-dependent four- and three-particle couplings, and F np 
and F npq are Q-dependent functions of the m 2 , depending on the particle types. The two-loop 
functions can be evaluated analytically using various methods developed in [15|-[19|. In general, 
one can write the results in terms of 10 basis functions, corresponding to the one-particle- 
irreducible connected vacuum graphs shown in Figure 1. These functions were given explicitly 
in ref. 11131 . 
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Figure 1: The one-particle-irreducible connected Feynman diagrams contributing to the two-loop ef- 
fective potential. Dashed lines denote real scalars, solid lines denote Weyl fermions carrying helicity 
along the arrow direction, wavy lines are for vector bosons, and dotted lines are for ghosts. The large 
dots between opposing arrows on the fermion lines in the FFS and FFV diagrams denote mass inser- 
tions. The FFS diagram is accompanied by its complex conjugate (the same diagram with all arrows 
reversed). The effects of the VV, VVV, and ggV diagrams can always be combined into a "gauge" 
contribution. The loop integral functions associated with these Feynman diagrams are given explicitly 
in ref. [13]. 
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In the MSSM, there are two Higgs VEVs, v u and v d . The evaluation of the two-loop ef- 
fective potential as a function of the VEVs reduces to determination of the relevant tree-level 
field-dependent couplings and masses. I will do this in the approximation of no Yukawa or 
flavor- violating couplings for the first two families of (s)quarks and (s)leptons. However, all 
remaining complex phases which cannot be rotated away are maintained. The resulting expres- 
sions presented here are complicated, but suitable for direct evaluation by computer programs. 

The rest of this paper is organized as follows. In section Q I describe necessary conventions 
and define some coefficients used in the calculation. Section ^ contains the expressions for the 
effective potential in the MSSM up to two-loop order. Section f| discusses supersymmetric limits 
of the effective potential. In section [|, I present two-loop results for the Higgs scalar anomalous 
dimensions and the beta function of the vacuum energy. These are necessary for checking the 
scale-independence of the effective potential. Section |6| treats a numerical example. 

2 Conventions and setup 

In this section, I list the necessary facts and conventions used in this paper. The MSSM is a 
softly-broken supersymmetric theory, with SU(3) C x SU{2) L x U(l)y gauge couplings g 3 , g, and 
g' . The last is related to the GUT normalized coupling gi by g' = ^3/5 g\. The superpotential 
of the MSSM is: 

W = uY u QH u - dY d QH d - eY e LH d + fiH u H d} (2.1) 

where H u and H d are the Higgs chiral superfields, and the Q, L are the chiral superfields con- 
taining left-handed quarks and leptons and u, d, e are those containing the conjugates of the 
right-handed quarks and leptons. The quark and lepton superfields carry a suppressed family 
index, so that Y u , Y d , Y e are 3x3 matrices in family space. Gauge indices for SU{2) L and 
SU(3) C are also suppressed, as in ref. ||. The Higgs mass parameter /i can have an arbitrary 
phase. The soft supersymmetry-breaking part of the Lagrangian is: 

-Aoft = (-MtBB + ^MzWW + ^Magg 

+ua u QH u - ddL d QH d - ea e LH d + bH u H d ^ + c.c. 
+m 2 H jH u \ 2 + m 2 H jH d \ 2 

+Q ] m 2 Q Q + L^m 2 L L + um 2 u u ] + dm 2 d d ] + em 2 e e ] . (2.2) 

where B, W, and g are the bino, wino, and gluino fields, and the same symbols are used for 
scalar fields as for the corresponding superfields. Here, a u , scalar cubic couplings in 

the form of 3 x 3 matrices in family space. The Higgs fields have soft supersymmetry-breaking 
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squared-mass running parameters rnf {u , m2 H d i an d b. The first two of these are necessarily real, 
and by convention b is taken to be real at the renormalization scale Q at which the effective 
potential is to be minimized. This can always be achieved by a suitable field rephasing, and 
ensures that the VEVs obtained by minimizing the full effective potential will also be real. The 
squarks and sleptons have running soft supersymmetry-breaking squared masses nig, m 2 L , m 2 , 
m^, and m 2 , which are 3x3 Hermitian matrices in flavor space. 

The gauge-eigenstate complex scalar doublet Higgs fields that are components of left-handed 
chiral superfields are called H u = (H+,H®) and H d = (H%,Hj). The electrically-neutral com- 
ponents have real positive VEVs v u and v d respectively. The field-dependent tree-level squared 
masses for vector bosons are then given by 

m 2 w = ^(v 2 u + v 2 d ); (2.3) 



The tree-level Higgs potential is 

v = A + (H 2 + <J(I^T + l^ + r) + (l/i| 2 + <J(I^T + l^ 



.2 



+ b(H:H d - H»H») +c.c. 

+\(9 2 + 9' 2 )(\H° U \ 2 + |#+| 2 - \H°\ 2 - \H d \ 2 ) 2 + 9 -\H:H°/ + H° U H?\*. (2.5) 
Here A is a running field-independent vacuum energy, which must be included to maintain 



renormalization-scale independence of the effective potential p0[-[23|. The neutral Higgs scalar 
tree-level squared masses are obtained by diagonalizing the matrices: 

(\^ + m 2 Hu + ^(3v 2 u -vl) -b - (g 2 + g' 2 )v u v d /2 \ 
* V ~b ~ (9 2 + g' 2 )v u v d /2 |/i| 2 + m 2 Hd + ^l(3v 2 - v 2 u ) ) ' 1 ' 1 

m\ - (W + ^ + ^M-i® b ) (2 7) 



'i 



which are written in the (Re[if°], Re[if °]) and (Im[iT°], Im[if^]) bases, respectively. The complex 
charge ±1 Higgs scalar tree-level squared masses are obtained by diagonalizing the matrix 



_2 I „/2 r, „2 



" = ~ l b + g 2 v u v d /2 \^ + m 2 Hd + ^vj + ^vl ] - U< 



which is written in the (H£, H d *) basis. The gauge-eigenstate fields can be expressed in terms 
of the tree-level squared-mass eigenstate fields as: 

m) = (S) + vWtfO + v^ o (^) (2 - 9) 
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and 

#f*)=^ ± (tf+)> ( 2 - 10 ) 

where G° and G ± are Nambu-Goldstone fields, h°, H°, A , and H ± are the physical Higgs tree- 
level mass-eigenstate fields, and 



it! a = Q a ; ifc = * Po ; = p± p± (2.11) 
are orthogonal matrices determined by the requirements that: 

K'm k R a = (f J 2 J; (2.12) 

flamioB ft = (™f ° ); (2.13) 



m£ 


° ) 







m 2 G0 


°. ) 





m\ J 













Rslml,R P± = "T J • (2-14) 



These equations define the tree-level squared mass eigenvalues for the Higgs scalar sector and 
the mixing angles a, fio and (5±. Here I use a notation in which c and s with a subscript indicate 
the cosine and sine of the indicated angle. So, in accord with the orthogonality of the rotation 
matrices, c a means cos a, etc. Similarly, in the following, C2 a is taken to mean cos(2a), and so 
on. Also, for future convenience, define coefficients 

k u h° = k dH o = c a ; k uH o = —k dh o = s a ; (2-15) 

Kg° = k dA ° = isp Q \ k uA o = -k dG o = ic^; (2.16) 

Kg+ = k dH + = sp ± ; k uH + = -k dG + = c p± . (2.17) 

Conventionally, c a , cp , and Cp ± are taken to be positive. Because the minimum of the effective 
potential is not a minimum of the tree-level potential, the angles (3q and (3± for the rotations in the 
pseudo-scalar and charged Higgs sector are distinct from each other, and from (3 = tan -1 ^/^) 
at the minimum of the effective potential. Some care should be used to distinguish these. Also, 
note that unlike the case in the ordinary Standard Model, m G o ^ m G ± at tree level. Note also 
that even with arbitrary CP- violating phases, there is no mixing in the tree-level squared masses 
of the neutral scalar (h°, H°) and the pseudo-scalar (G°, A ) sectors, because of the freedom to 
choose the argument of b to be at any particular renormalization scale Q. 
In section |3|, V] and will appear, denoting sums over the lists 

<t>o <j>± 

0° = (h°, H°, G°, A ) and ± = (G ± , (2.18) 

respectively. 
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The neutralinos (A^; i = 1,2,3,4) and charginos (Cf 



1, 2) are mixtures of the elec- 



troweak gaugino and Higgsino fields. In the {B, W°, H%, H®) gauge-eigenstate basis, the neu- 
tralino mass matrix is 



( M x -g'v d /V2 g'v u /V2\ 

M 2 gv d /y/2 -gv u /V2 

-g'v d /V2 gv d /V2 -fi 

V g'v u /V2 -gv u /V2 



(2.19) 



This is diagonalized by a unitary matrix N: 

N*Mf f N~ 1 = diag^ m^ m^ ), 



(2.20) 



where the mass eigenvalues rrif f _ are all real and positive. This can always be accomplished (for 
any phases of Mi, M 2 , and //) by the following procedure. First, define E to be a matrix whose 
columns are orthonormal eigenvectors of the Hermitian squared-mass matrix M^Mfj, arranged 
in order of increasing eigenvalue. The orthonormality means that E^E = 1. It follows that 

(2.21) 



E T MffE 



P 2 m D , 



where is the matrix on the right-hand side of eq. ( |2.20| ), and P is a diagonal phase matrix. 
Then 



N = PE^ 



[2.22] 



satisfies eq. ( |2.20|) , and also 

NM^M^N- 1 = diagl 



m% x ,ml 2 ,m\,m\ A ). 



;2.23) 



This procedure gives the tree-level field-dependent neutralino squared masses and mixing matrix 
Nij. (A suitable generalization of this method can be used to transform the mass and squared- 
mass matrices into a real positive diagonal form for any fermions in any theory.) 
The tree-level chargino mass matrix is given by 



Mc 



(M 2 gv u 

\gv d /i 



U*M e V* 



and is similarly diagonalized by unitary matrices U and V according to: 

mA 
m c 2 

where again mg. are real and positive. Then 

' m% 



VMtMcV- 1 



UM^M^U- 1 



Ci 







(2.24) 



(2.25) 



(2.26) 



The tree- level squared mass of the gluino, denoted m~ = |M 3 | 2 , does not depend on the 
VEVs v u and Vd- 

In most of this paper, I will employ the approximation that only the third-family Yukawa 
couplings y t , yb, and y T are significant, so that: 



Y 



u 





















(0 







(i 





'')■ 


Y d = 










Y e = 








: ) 







yt) 




\0 





VbJ 




\0 





Vr) 



(2.27) 



By suitable field redefinitions, they can be chosen real and positive. The non-zero quark and 
lepton squared masses are therefore: 

222 222 222 /<-> <->o\ 

m t = y t v u , m b = y h v d ] m T = y T v d . (2.28) 

(Effects of the other Yukawa couplings are certainly smaller than the dominant 3-loop order 
contributions.) 

I will also assume here that the soft supersymmetry-breaking scalar cubic couplings involving 
first- and second-family sfermions vanish, so 






ad = ; a e = , (2.29) 



and that, in the same basis, the soft supersymmetry-breaking scalar squared mass parameters 
are diagonal in family space, so 

lm\ x \ 
m 2 Q = m% 2 , etc. (2.30) 
V m\J 

The parameters^ at, ^b, and a T can in general be complex, and 2 3 etc. are real soft running 
OR' squared mass parameters. The squared masses of the first-family squarks and sleptons are 
then given by: 



m 



k = m Q, + A ^ m \ L = m Q, + A d L ; (2-31) 

< R = < + An fl ; m\ R = m 2 di + A d - K ; (2.32) 

ml = m 2 Li +A, e ; m\ = m\ x + A Sl ; (2.33) 

m\ R = m 2 ei + A iR , (2.34) 



The D-term contributions are: 



1 



A/ = l{If9 2 -Y } g a ){vl-vl), (2.35) 

with I j and Yj defined to be the third component of weak isospin and the weak hypercharge of 
the left-handed chiral superfield containing the squark or slepton /: 

^In the literature, one often sees rescaled quantities A t — at/Ut> — db/yt, A T = a T /y T , and B = bj \i. 



9 



1/2 -1/2 1/2 -1/2 
1/6 1/6 -1/2 -1/2 -2/3 1/3 1 



Exactly analogous expressions hold for the second-family squark and slepton squared masses 
m? L , m 2 L , rn\ R , m 2 R , m 2 , rnfi L , and nv~ H , with the subscripts 1 replaced by 2 on the right-hand 
sides of eqs. (^31|)-(^). 



For the stop, sbottom, stau, and stau-neutrino squared masses: 

m ? = ( m Q 3 + Vt v l + A u L v u< ~ v d /jy t \ . , 2 3g x 

* ' \ v u a t - v d n*y t m 2 U3 + y 2 v 2 u + A^ ) ' 

m ? = ( ™ 2 q 3 + vW d + Ai L v d a* b - v u [iy h \ . ^ 3? * 

b \ v d a b - v u fi*y b m 2 d3 + y%vj + A dR J ' 

f " " V - v u n*y T m 2 e3 + y 2 vj + A e - B J ' 

«l = m| 3 + A, e . (2.39) 

Here the stop, sbottom and stau squared mass matrices are given in the /#) bases. The 
eigenvalues of eqs. ( |2.36| )-( |2T3"%| ) give the squared mass eigenvalues m~ , m| , and m? for i — 1,2. 
The squared-mass matrices are diagonalized by unitary transformations: 

for f = t,b, r, with 



so that 

V^.-=(1- m \). (2-42) 

and similarly for the sbottoms and staus. Unitarity of the matrix Xj allows one to write 
Lf = R* ; = Cf, and Rf = —L% = s?, with 

| C/ ~| 2 + | S/ -| 2 = 1. (2.43) 

(If the off-diagonal elements of the squared mass matrix are real, then cj and sj are the sine 
and cosine of a sfermion mixing angle.) In the following, I will also make use of coefficients: 

xj = Ij, (/ = first or second family sfermion); (2.44) 

T^l-^l 2 ' x li = ~ o^J 2 ' = 2' x h = ~-j\Lfi\ 2 1 (2-45) 



x 



u 
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and 



x'j = Yj, (/ = first or second family sfermion); (2.46) 

X L = h L . |2 _ l\ R |2 xL = I| r |2 + I )jR |2 (2>47) 

< = -^l^l 2 + l^l 2 ; < = -^- ( 2 - 48 ) 

When V] appears in section it will denote a sum over the list 
/ 

f = (ul, d L , u R , d R , u e , e L , e R , c L , s L , c R , s R , u^, p, L , J1 R , t u t 2 ,bi,b 2 , u T , n, r 2 ). (2.49) 
The symbol 

n f -={ h l = Slept ° n (2.50) 
[3, / = squark 

denotes the number of colors. 

To summarize, evaluation of the effective potential can proceed as follows. At a renor- 
malization scale Q, choose values for the set of input data consisting of the 33 DR' running 
parameters: 

v u , v d , (2.51) 

93, 9, 9', Vt, Vb, Vt, (2-52) 

m Q t ' m h m li' m % m lv (» = X ' 2 ' 3 ) ( 2 - 53 ) 

m 2 H u , m H d , b , M, (2-54) 

M 3 , M 2 , Mi, at, a b , a T , (2.55) 

of which the last 7 may be complex. Using the preceding protocols, evaluate the mixing param- 
eters 

a, fa, P ± , (2.56) 

Nij, U tj , % (2.57) 

t>t Rti, L~ bl , Rli> L Ti, Rfi, (2.58) 

and the 39 distinct field-dependent tree-level squared masses: 

m 2 w , m|, (2.59) 

2 2 2 2 2 2 /n rn\ 

m h o, m H o, m G o, m A o, m G ±, m H ±, (2.60) 

m % m \s m cv ( 2 - 61 ) 
11 



m t ,m b ,m T , (2.62) 

2222222 /n ri\ 

ul ' cZl ' «fl ' d R ' ' e L i e R > V / 

2222222 /<-> 
2222222 /<-> r>r\ 

m- tl , m h , m Si , m^, m^, m fl , m f2 . (2.65) 

These squared masses (and implicitly Q) then become arguments for the functions Fsss, Fss, 
Fffs, Fws> F ssv, F vs , F vvs , F ffv , F-p^y, and F gauge appearing in the next section, and 
defined explicitly in eqs. (2.17)-(2.22), (4.12)-(4.16), and (6.21)-(6.30) of ref. 0. (For the sake 
of uniformity of notation, I use the symbols Fsss, Fss, Fffs, Fpp S , an d Fssv in place of fsss, 
fss, fFFS, fjTs' an d fssv for the functions which are the same in the MS and DR' schemes. 
All two- loop functions used in the present paper are DR' ones.) To simplify the notation, I 
adopt the convention that the name of a particle is synonymous with its squared mass when 
appearing as an argument of one of these functions. So, for example, F S ss(h°, A , G°) means 
F S ss{m 2 h o,m 2 AO ,m 2 Go ). 

3 MSSM effective potential 

3.1 Tree-level and one-loop contributions 

The tree-level contribution to the effective potential for v u and Vd is: 

V<°> = A + (H 2 + Ol + + miy d - 2bv ^ d + \{g 2 + g' 2 )(vl - v 2 f. (3.1) 
The one-loop contribution in the DR' scheme is: 

V® = ^M0°) + 2^/ i (0 ± ) + 2^n / / i (/)-2^/ i (iV,)-4^Ma) 

cj>o 4>± J i=l i=l 

-I6h(g) - I2h(t) - 12h(b) - 4/i(r) + 3h{Z) + 6h(W), (3.2) 

where h(x) is the function defined in eq. ( |L3| ), and the name of each particle is used to denote 
its squared mass. 

3.2 (S'S'S'-diagram contributions 

In this subsection, I list the contributions to the MSSM two-loop effective potential from dia- 
grams with three scalar propagators. These all involve the function Fsss{%, V, z ), with arguments 
x, y, z equal to tree-level field-dependent scalar squared masses. 
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The contributions from diagrams with three 

V d>°6°6 — 



(g 2 + g' 2 ) 2 

32 



Higgs scalar propagators are: 
(c a v u + SaVaf^Fsssih , h°, h°) + 2s 2 2(3o F S ss(h°, A , G°) 
+4p [Fsss(h°, G°, G°) + F S ss(h°, A , 

+(s a v u - c a v d ) 2 ^cl a F sss {H^ H°, H°) + 2s^ F sss (H°, A , G°) 

+c% [F S ss(H°, G°, G°) + Fsss(H°, A , 
+[s a (l - Qc 2 a )v u + ca(l- Qsl)v d ] 2 F S ss(h°, h°, H°) 
+ [s a (l - 6c 2 a )v d - c a (l - 6s 2 a )v u } 2 F sss (h , H°, H°) 



(3.3) 



and 



V ( 



(2) 



1 



g 2 c 2 /3 ± {s a v u - c a v d ) + g' 2 s 2 j3 ± (c a v u + s a v d ) F SS s(h°, G ± , F ± ) 



+2 
+ 



g' 2 C2p± (c a v u + s a v d ) - g' 2 s 2f3± (s a v u - c a v d ) F S ss(H°, G ± , H 



g 2 {c a v u - s a v d ) + g 2 s 2l 3 ± {s a v u - c a v d ) - g' 2 c 2l3± {c a v u + s a v d )\ F sss (h , G ± , G ±X; 

i 2 

g 2 (c a v u - s a v d ) - g 2 s 2/3± (s a v u - c a v d ) + g' 2 c 2l3± (c a v u + s a v d ) F SS s(h°, H ± , H 



+ 
+ 
+ 



g 2 {s a v u + c a v d ) - g 2 s 2/3± (c a v u + s a v d ) - g' 2 c 2f3± (s a v u - c a v d )\~ F SS s(H°, G ± , G ± ' 
g 2 (s a v u + c a v d ) + g 2 s 2f3± (c a v u + s a v d ) + g' 2 c w± (s a v u - c a v d ) F SSS (H°, H ± , H 



+2g\c Po v u - s^fFsssiG , G ± , ±) + 2g\sp Q v u + c Po v d ) 2 F sss {A^ G ± , H*) \. (3.4) 



o n± 



The contributions from diagrams with two sfermions and a neutral Higgs scalar are: 



V 



(2) 



0°//' 



E \ 2F sss(<f>°, f,f). 



(3.5) 



0°,/,/' 



The couplings involving first- and second-family sfermions, and the tau sneutrino, are non-zero 
only when / = /'. They are given by: 



V/7* 



-y= (ijg 2 - Yjg' 2 ^j Re[k U(j ,ov u - k d(/> ov d ], 



(3.6) 



with Ip Yj defined in section 0. The remaining couplings for the third-family sfermions are: 



A 



<t>°ut* 



L iiLiM°u L ui + R u R tK u R u R - ^v u y 2 Re[k u ^]{Li.L*z + R^Rl 



'u L u L i J - l/ t i J -"t :j "(P"u R u R 



--j={Kp<h - k d ^^*yt)L k Rl - -j={k* u ^a* t - k^ofxy^R^. 



(3.7) 



V2v d y 2 b MW]( L lLt +R- b Rt 
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-^{k d4>0 a b - k* w pH*y h )L i .R* j - — (fc^aj - k u4> ^y b )R~ h L\-, (3.8) 

1 1 

-^(W' a ^ - kl jfl n*y T )Lf i R^. - -j=(k* d4>0 a* T - k u<t> ^y T )R fi L* f] . (3.9) 

The contributions from diagrams with two sfermions and a charged Higgs scalar are: 

V$ rr = E n^j^Fsssi^JJ'). (3.10) 
4>+JJ' 

Here the non-zero couplings involving the first- and second-family sfermions are: 

g 2 

^<j>+e L v* = ^cj>+ii L v* = \+d L u* L = ^<P+s L c* L = —(k U( p+V u + k d( p+V d ), (3-11) 

and those involving the third-family sfermions are: 

Vm* = Vl { x <t>+d L u* L + y 2 t v uK<t>+ + y 2 b v d k dlj)+ ) + Ri i B%.y t y b {k dll> +v v + k u4> +v d ) 

+L~ b .R~. (k u4> +at + k d< p+/j,*y t ) + R b .L~. (k d<t>+ a* b + k u4> +/j,y b ); (3.12) 
V-fii>; = L fi (^+e L i>* + yr v dkd<i>+) + RfAk^+a* + k U(j} +iiy T ). (3.13) 

3.3 ^-diagram contributions 

In this subsection, I list the contributions to the MSSM two-loop effective potential which come 
from diagrams with two scalar field propagators. They all involve the function F S s(x,y). 
The contributions proportional to g\ are: 

{2 
E [\ L uL% ~ RkRi/FssiUjj) + \L b L\ 3 - \RlfFssiKbj) 

+F ss (u L , u L ) + F ss (ur, u r ) + F ss (d L , d L ) + F ss (d R , d R ) 

+F ss (c L , c L ) + Fss{c R , c R ) + F ss (s L , s L ) + F ss (s R , s R ) J. (3.14) 

The contributions from diagrams with two Higgs scalar propagators are: 

V$lo = g -^f{3cUFss(h°, h°) + F SS (H°, H°)} + 3ci, [F ss (A° , A") + F SS (G°, G )] 
+(4 - 6cl a )F ss (h , H°) + (4 - ^)F SS (A°, G°) 

+2c 2a c 2Po [Fss(h ,A )+F S s(H ,G )-F ss (h ,G )-Fss(H ,A )]y, (3.15) 
V^ ± = ^^{c^jF^tf^+F^^ (3.16) 
V $fi = l^^-^^+g^c^Fssih^H^ + FssiH^G^] 
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+ [g 2 (l + s 2a s w± ) - g' 2 c 2a c 2h }[Fss(h , G ± ) + F SS (H°, 

+ [g 2 (l - s 2A) s 2/3± ) + g' 2 c Wo c w± ] [F SS (A°, if ±) + F SS (G°, G ± )] 

+[<? 2 (1 + s Wo s w± ) - g' 2 c Wi) c w± ] [F SS (A°, G ± ) + F SS (G°, ii ± )]}. (3.17) 

The contributions from diagrams involving electroweak gauge couplings and one sfermion 
and one Higgs scalar propagator are: 



v?/*' = 7EEn/(V-^)(IWI 2 -l^l 2 )^5(/,/); 



(2),flV 2 



E E »/( V + *> /2 ) - 



(3.18) 
(3.19) 



The contributions from diagrams involving electroweak gauge couplings and two sfermion 
propagators are: 

,2 



V K ~~ 
ff 



{2W = y \ E npf.XfXf.Fssd f) + E n F )F ss {f\ /) 



/,/' 



/ 



+ ol^i^r^s(*i,* 2 ) + -I^LgJ^^^s) + -|L fl L f2 | 2 F 5S (n,f 2 ) 

2 



2 
+3 



Fss(u L ,d L ) + F ss (5l,S l ) + E l L *> L L | 2 *ss(*i, &j) 



2 _ I 

+^S5(eL, £ e ) + F ss (fi L , Up) + ^ |L f J 2 F55(ri, z/ T ) h 

1=1 J 



(3.20) 



T/(2),9' : 
// 



,/2 



E n f ~n fl x' f ~x' fl Fss(f, /') + E n f ~x]F ss (f, f) 



l 



+ -|i^? 2 - 4^^!^^,^) + -\L~ b Ll + 2^^^55(61,62) 



6 

+i|L fl L; 2 -2 J R flJ R; 2 | 2 F55(r 1 ,f 2 ) 



(3.21) 



The contributions from diagrams involving Yukawa couplings and one sfermion and one Higgs 
scalar propagator are: 



V 



m,y 2 



f<t>° 



EE{|l/* 2 IM 2j? ss(*i,/) + ItiMFssM ) + ^\k d ^F ss (nA )}; (3.22) 



vfX = E[E{% 2 K~J 2fc V + v?| + y 2 r\Rn\ 2 k 2 d<t>± F S s(n,<P d 



i>± L i=l 



+3(y?\L ii \X (l>+ +y 2 b \R l fk 2 d(t>+ )F ss (b l , ( p ± )} + y 2 T k 2 dlP± F ss (V T , (3.23) 
The contributions from diagrams involving Yukawa couplings and two sfermion propagators are: 

T/(2),2/ 2 



E y'MLifih I 2 + ^RlRf^FssiU, tj) 

i,j=l 1 
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+y 2 b (S\L ii R ij | 2 + QL^RlR^Fssih, b 3 ) + V ^\L fi R fj + R fi L f] \ 2 Fss(n, fj) 
+3(y 2 \R k L i . | 2 + y 2 b \L h R h | 2 )F 55 (^) + G^Re^i?? R, .L^Fssih fj) 

+y 2 T j2\RffF ss (T l ,v T ). (3.24) 
i=i 



3.4 FFiS-diagram and FFS'-diagram contributions 

In this subsection I list the contributions to the MSSM two-loop effective potential from diagrams 
involving scalars and fermions. These include the FFS diagrams (without chirality-ffipping 



fermion mass insertions) and the FFS diagrams (which do have two such fermion mass insertions 
on different propagators). They therefore involve the functions F FF s(x,y, z) and Fp-p S (x,y, z). 
The contributions from diagrams involving the gluino are given by: 

r 2 

V m = 8 ^3 E [ F FF,s(t, g, U) - 2Re[L ii Rl]m t m~ g F- FFs (t, g, U) 

+F FFS (b,g,bi) - 2Re[L i .Rl i ]m b rrigF ¥Fs (b,g,b i ) 

+F FFS (0, g, ml) + F FFS (0, g, d L ) + F FFS (0, g, u R ) + F FFS (0, g, d R ) 

+F FFS {0,g,c L ) + F FFS {0,g,s L ) + F FFS {0,g,c R ) + F FF5 (0, ^, s R )|. (3.25) 

The contributions from diagrams involving neutralinos, chiral fermions, and sfermions are: 
4 f 2 r 

Vf&f = E E 3(|F t ^| 2 + |y w /)^^^ 

i=i {j=i 3 

+3(\Y b Nr b *] 2 + \Y mii f)F FFS (b,^ 

HKn^ I 2 + Kn^ \ 2 )Fffs(t, Ni, fj) + 2Re[Y T ^Y-^ Jm^F^lr, N u r,) 

3 J 3 J 

v e ) + F FFS {0, Ni, z/ M ) + F FFS {0, Ni, v T )\ 
+3\Y uffiill \ 2 [F FF s(0, N h u L ) + F FFS (0, N t , c L )} 



+31^^11^5(0, Ni, u R ) + F FFS (0, Ni, Z R )\ 
+2>\Y dfl ^ \ 2 [F FFS (0, N h d L ) + F FFS (0, N u s L )} 
+3\ Y dNj R \ 2 l F FFs(0, N h d R ) + F FFS (0, N h s R )} 
MYeN.ei \ 2 [Fffs(0, N t , e L ) + F FFS (0, N h Ji L )} 

+ \Y- e ^ eR \ 2 [F FFS {U,N l ,e R ) + F FF s(0, Ni, /iR)]j. (3.26) 
Here, the fermion-neutralino-sfermion couplings are: 

y u N iK = -^9N* 2 + ^N* ± y, Y^- r = ^g'N*,; (3.27) 
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*W L = 7f (s^ - f^A); *W« = -jM; (3-28) 

n^e-i = + *Wi); n*^ = -^V^A; ( 3 -29) 

W = -^(-<?^ + </^A); (3-30) 

Y tNit* = L iYuNiU* L - R % N iAVt\ Y tNdj = RijYuNiUR ~ ^t^tiVu (3.31) 

Ymr, = L l Y d Ndl ~ %*aW, Y mh = R~ b Y mjR - L b N* 3 y b ; (3.32) 

Y Tf!if; = L%Y e ^ L - R* f .N* 3 y T ; = ll : ) TX , !: - L f] N* 3 y T . (3.33) 

The contributions from diagrams involving charginos, chiral fermions and sfermions are: 



i=l kj=l 



+3(lncV* I + \ 2 )F F Fs(b, Q, tj) + 6Re[y b< 5 if *y^ t -.]m 6 m (5i F rF5 (fo, Q, t 



j j 



+H\Y tC w I 2 + \ Y m \ 2 ) F ^s(t, Q, bj) + QReiY^^r jmtmc^sit, Q, b d 

3 J 3 J 

HKc^l 2 + Kc^\ 2 )Fffs(t, Ci, V T ) + 2Re[Y Tdi0t Y Tdi ~ T ]m T mc i F ms (T, Q, V T ) 
+3|y^ ifi . \ 2 [F FFS (0, Ci, u L ) + F FFS (0, Ci, c L )\ 
+3\Y uCii \ 2 [F FFS (0, Q, d L ) + F FFS (0, Q, s L )} 



+\Y e c i0i \ 2 [Fffs(0, C i} v e ) + F FFS (0, C h u, 



+\Y l/edi , l \ 2 [F FFS (0,C i ,e L ) +F FFS (0,Ci,ji L )]j. (3.34) 
Here, the non-zero fermion-chargino-sfermion couplings are: 

YdCiU* L = Y eC t u* = Y rCiU* = ~9V*1' Y uCid* L = Y v e C t e* L = ~9^*\\ (3.35) 

Y^ = -L* ijg V* + R\y; 2 y t ; Y m = L h U* 2 y b ; (3.36) 

Y tCib* = -fyUti + RtU*2y b ; Y m , = L bj V*y t ; (3.37) 

Y VT&in = -L^gUZ + R* fj U* 2 y T ; Y^ = y T U* 2 . (3.38) 

The contributions from diagrams involving Higgs scalars and Standard Model fermions are: 

V f%o = ^Ei 3 ^ 2 [l^r^F5(t,t,0°) + (A;^o) 2 m t 2 F Ws (t,t,0 )] 
1 4,0 L 

+Zy 2 b [\k d<t>0 \ 2 F FFS {b, b, 0°) + {k^fmlF^ib, b, /)] 

+y 2 l\k d(f)0 \ 2 F FFS (T,T,<j ) ) + (A; # o) 2 m 2 F Ws (r,r,0°)]|; (3.39) 
^//'V = J2) Hyt k l<t>+ + yl k %+) F FFs(t, b, ± ) + Qk u4> +k d<t) +y t y h m t m b Fpp S (t, b, ± ) 



0± 

+y 2 A + F FFS ^,r,^)\. (3.40) 
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The contributions from diagrams which involve a Higgs scalar propagator and two chargino 
and/or neutralino propagators are: 

vgL = EEll^c-^r^5(Q,Q,0°) 

+Re[Y e+d -^Y d+d -^}m e m d Fj^siCi, C h 0°) }; (3.41) 
+Re[(Y m ^) 2 }m ff m^F- FFs (N l , Nj, 0°)}; (3.42) 

i=i j=i 4>± k 

+2Re[y^ J v.^y 6 - J v.^ + ]m (5i m J v.Fpp 5 (C i , % 0±)}, (3.43) 



where the necessary couplings are: 

Ycfc-o* = -^[k^V^ + K^U*,]-, (3.44) 

*W = ^(^-^a)(^oA7 4 -^oA7 3 ) + (^i); (3-45) 
^iv,0- = K < p + [9V* 1 N* 4 + -Lv^gN^ + g'N*,)]; (3.46) 

= W[gU; i N; 3 -^=U* 2 (gN; 2 + g'N* 1 )]. (3.47) 
3.5 iS^y-diagram contributions 

In this subsection I list the contributions to the MSSM two-loop effective potential coming 
from diagrams involving one vector and two scalar propagators. These all involve the function 

Fssv(x,y,z). 

The contributions from diagrams involving the gluon and squarks are: 

V Sl = 2glJ2Fssv(q,q,0). (3.48) 

The contributions from diagrams involving the photon and sfermions are: 

V}% = 2{g (f g a ) ^ n fQ 2 f F ssv(llO)- (3.49) 



Here Qj denotes the electric charge of the sfermion. 
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The contributions from diagrams involving the W, Z bosons and sfermions are: 
1 



(3g* - g y[Fssv(u L , u L , Z) + F ssv {c L , c L , Z)\ 



ffz 24(g 2 + g' 2 ) k 

+(3/ + g' 2 ) 2 [F SS v(d L J L , Z) + F S sv(s L , s L , Z)\ 

+16g' 4 [F S sv(u R ,u R , Z) + F ssv (c R , c R , Z)\ 

+4g' 4 {F ssv (d R , d R , Z) + F ssv (s R , s R , Z)\ 

+3(<? 2 ~ 9 ,2 ) 2 [Fssv(e L , e L , Z) + F SSV (fa, fa, Z)\ 

+3(# 2 + g' 2 ) 2 {F S sv(v e , u e , Z) + Fssvipp, Vp, Z) + F SS v(vt, v T , Z)\ 

+l2g ,4 [F ssv (e R , e R , Z) + F ssv (fa, fa, Z)\ 

+ W-9 l2 )L i M j -^' 2 ^M j \ 2F ssv{u,t j ,z) 

+ £ W + 9' 2 )L~ bi Ll -2g' 2 R ii Rl \ 2 F SS v(b i ,b j ,Z) 

+3 £ \{9 2 ~ g' 2 )L h L% 3 -2g' 2 R ft R; j \ 2 F ssv (T u T J ,Z)\; 
i,j=i ) 



(3.50) 



V-. 



(2) 



2 r 2 

3F 55y (5 L , ck, WO + 3F 55y (5 L , sl, + 3 £ |L t - L s . \ 2 F SSV {U, bj, W) 

2 



+F ssv (u e , e L , W) + F S5y (^, fa, W)+Y. \Lr t \ 2 Fssv(vT, %, W) 



i=i 



(3.51) 



The contributions from diagrams involving Higgs scalars and electroweak gauge bosons are: 



1/(2) 

v </>±</>±7 ~~ 


g 2 g a 


2(^ + g?) 


T/(2) 


{g 2 -g' 2 ) 2 
Kg 2 + g' 2 ) 


T/(2) 


g 2 + g' 2 \( 
8 1 (< 



FssviH*, H ± , 0) + F SSV (G ± , G ± , 0) 
Fssv 

(H ± ,H ± ,Z) + F SSV (G ± ,G ± ,Z) 
c aC/3o + s aS/3o ) 2 \F ssv (h°, A , Z) + F SSV (H°, G°, Z) 



v ,j) (f, ± W 



+(s aC/3o - c a s Po Y [F ssv (h°, G°, Z) + F SSV (H\ A\ Z) 

j{(c aC/3± + s aSf5± ) 2 [F S sv(h°, H ± , W) + F SSV (H°, G ± , W) 

+ {s a c p± - c a s p± ) 2 [F ssv (h°, G ± , W) + F SSV (H°, H ± , W) 
+(c A)C/3± + Sf3o s p± ) 2 [f ssv (A°, H ± , W) + F SSV (G°, G ± , W) 

+(s A)C/3± - c A)S/3± ) 2 \F SSV (A°, G ± , W) + F SSV (G°, H ± , W) 



(3.52) 
(3.53) 

(3.54) 



(3.55) 
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3.6 V^-diagram contributions 



In this subsection, I list the contributions coming from diagrams with one vector and one scalar 
propagator. These all involve the function F vs (x,z). This function vanishes when the vector 
squared mass variable x is zero, so only the W and Z bosons contribute. 

The contributions from diagrams with one electroweak gauge boson and one Higgs scalar 
propagator are: 



(2) 



W4> 



£ £ F VS (W, 0°) + ^ E Fvs(W, ± ); 



V?l = 9 -^^F V s(Z,<f>°) + 



z<t> 



8 



(g 2 - g' 2 ) 2 
±{g 2 + g' 2 ) 



(3.56) 
(3.57) 



The contributions from diagrams with one electroweak gauge boson and one sfermion prop- 
agator are: 



V {2 1 = 

Wf 



v {2 2 = 

Zf 



sf ty En f \x f \F vs (Wj) ] 
f 

1 



(3.58) 



16«/ 4 F vs (Z,u R ) + F vs (Z,c R ) + |i? t -J 2 F y5 (Z, t,) + \R i2 \ 2 F vs (Z,t 2 ) 



12(<? 2 + <?' 2 , , 

+ (3g 2 - g' 2 ? [Fvs(Z, u l ) + F VS (Z, c L ) + \L h \ 2 F VS (Z, h) + \L i2 \ 2 F vs (Z, t 2 ) 
+4g' i [F vs (Z,d R ) + F vs (Z,s R ) + \R~ hi \ 2 F VS (Z, h) + \R h \ 2 F VS (Z, b 2 ) 
+(3g 2 + g' 2 ) 2 {F vs (Z,d L ) + F vs (Z,s L ) + \L- bl \ 2 F VS (Z, h) + \L h \ 2 F vs (Z, b 2 ) 
+ 12g' A [F vs (Z,e R ) + F vs (Z,ji R ) + \R fl \ 2 F VS (Z, f ± ) + \R f2 \ 2 F vs (Z,r 2 )\ 
+3(<7 2 - g >2 ) 2 [F vs (Z,e L ) + F vs (Z,fi L ) + lUfFvsiZ,^) + \L f2 \ 2 F VS (Z, r 2 ) 

+3(g 2 + g' 2 ) 2 \F VS (Z, V e ) + F VS (Z, u„) + F VS (Z, V T ) 



(3.59) 



3.7 VV ^-diagram contributions 

The contributions from diagrams with one scalar and two vector propagators are: 



Vi 



(2) 



2„/2 



9 9 



{{c p± v d + S(i± v u ) 2 [g 2 F vvs (W, 0, G ± ) + g ,2 F vvs (W, Z, G ± ) 



2(g 2 + g' 2 ) 

+(cp ± v u - s p± v d ) 2 [g 2 F vvs (W, 0, if ±) + g' 2 F vvs (W, Z, if ±)] } 
+\{c a v u - s a v d ) 2 [{g 2 + g' 2 ) 2 F vvs (Z, Z, h°) + 2g 4 F vvs (W, W, h°) 
+ \{c a v d + s a v u ) 2 [{g 2 + g >2 ) 2 F vvs (Z, Z, ff°) + 2g 4 F vvs (W, W, H )} . (3.60) 
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3.8 FFF-diagram and FFF-diagram contributions 



In this section I list the contributions from diagrams with one vector boson and two fermion 
propagators. These include cases with no chirality-flipping mass insertion [involving the function 
Fffv(x, U, z)] and those with two such mass insertions on different propagators [involving the 
function F- ¥Fv (x, y, z)\. 

The contributions from diagrams with two Standard Model fermions and one vector boson 

are: 

V£l = 4 ^3 [F F Fv(t, t, 0) - m 2 t F m {t, t, 0) + F FFV (b, b, 0) - m 2 iW&, b, 0)] ; (3.61) 
r{2) ~ j^2^ F ^v{t^)-m 2 t F- WPv {t,tM 



V (Z> 

v ffi ~ 



V, 



(2) _ 



9 +y ^ 

+^[F FFV (b,b,0)-mlF JTv (b,b,0)}+F FFV (T,T,0)-mlF ¥Fv (r,T,0)^ (3.62) 
(51g 4 + Qg 2 g' 2 + 83g' 4 )F FFV (0, 0, Z) 



" z " 24(^ 2 + g' 2 ) L 

+ (9g 4 - Qg 2 g' 2 + 17g' 4 )F FFV (t, t, Z) + 8g' 2 (3g 2 - g^F^t, f , Z) 

+ (9g 4 + 6g 2 g' 2 + 5g' 4 )F FFV (b, b, Z) + 4g' 2 (3g 2 + ^ 2 )m 2 F PFy (6, 6, Z) 

+(3(? 4 - 6g 2 g' 2 + 15«/ 4 )F FFy (r, r, Z) + 12«/V - ^^^(t, r, Z) 



V, 



. (2) 9* 

//w - y 



3F FW (t, b, W) + F FF y(r, 0, W) + 8F FFy (0, 0, W) 



(3.63) 
(3.64) 



The contributions from diagrams involving charginos and/or neutralinos and electroweak 
vector bosons are: 

l2 E ^FFv{Ci, Ci, 0) — m~_Fpp V (Ci, Ci, 0) ; 
~r 9 ,:=i L 



V< 2 > 
NNZ 



ccz 



9 2 


+ g' 2 


9 2 


+ g' 2 




2 






9 2 


+ g' 2 


2 



E 

2 

E 

£.7=1 



0;; L | 2 F FFy (iV l? iVj, Z) + {O'^fm^F-^yiN^ N j: Z) 

:\o'g\ 2 + \o'*\ 2 )F F MCiA,z) 



(3.65) 
; (3.66) 



-20' t fO' t f*m 5 m 5 F mv (C t , C v Z) 

v^ w = ^EE{(lo§l 2 + |oSl 2 )W^c,-,wO 



(3.67) 



i=l j=l 

where the necessary couplings are: 
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(3.68) 



(3.69) 



o't = -y^-\^ + -r^r 2 ^ (3-70) 



o'if = -u* 1 u jl -lu* 2 u j2 + -^- 2 8 ij] (3.7i ; 

1 



OS = N i2 V; x -^=N^V; 2] (.3.72) 



Og = KUji + ^=KUi2. (3-73) 



There is also a contribution involving gluon and gluino propagators: 

Vff g = 12g 2 3 [F FFV (g,g,0) ~ |M 3 | 2 F^ y (^, ^, 0)] . (3.74) 

This is independent of the VEVs v u and Vd, and therefore is irrelevant to the minimization 
of the effective potential. However, it is Q-dependent, and therefore must be included when 
renormalization group consistency is checked. 

3.9 Pure gauge contributions 

The contributions involving only vector and ghost fields are 



V £L = o, 2 9 , m k 2i? g auge(^, W, 0) + y 2 F gauge (H/, W, Z)\ . (3.75) 



gauge 2 ^ 2 + gf2 ^ 

This concludes the list of contributions to the two-loop effective potential in the MSSM. 
Partial results for the two-loop contributions in the approximation that g, g', y T , and a T vanish 
and there are no CP- violating phases had previously been given in P|-[IT]. (Applications of 



these results to the Higgs scalar boson mass spectrum have been made in refs. |p4||-|p7||.) The 
two-loop contributions to the effective potential involving only Standard Model fields were given 
in Hl7| , but in the MS scheme, which is not convenient for the supersymmetric extension. 



4 Supersymmetric limits 



The results of section |3] can be checked by considering non-realistic limits in which supersymme- 
try is restored. Unbroken global supersymmetry requires that the effective potential vanishes. 

One such limit occurs if all supersymmetry breaking parameters are 0, and v u = = 0. The 
only massive particles in the theory are then the members of the Higgs supermultiplets, with a 
common squared mass x = \fi\ 2 . The one- loop effective potential is easily seen to vanish, since 
it is a supertrace over the squared masses. From the results of section |3], one finds that: 

V W = ^ p [F ss (x,x)+AF FFS (0 ,x,x) + F ssv (x, x, 0) + F FFV (x, x, 0) - xFpp V (x, x, 0) 
+2(3yf + 3y 2 b + y 2 T ) \xF SS s(0, 0, x) + F FFS (0, 0, x) + 2F FFS (0, x, 0)1 . (4.1 
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Each of the quantities in brackets indeed vanishes, using the expressions in ref. [|Tj| . 

Another way of maintaining supersymmetry is to again make all supersymmetry-breaking 
parameters 0, but now also require \i = 0, and take the VEVs along a .D-flat direction v u = 
Vd- The parameters g, g', gs,yt,yb,y T remain arbitrary. Then one finds that = = 0. 
However, the two-loop contribution from section [| does not vanish: 

2 



V (2) 



g 2 + g° 



3J{m;) - 3J(mt) - J{m l T ) 



(4.2) 



where 



J(x) 



x 



ln(x/Q 2 



(4.3) 



This might seem to violate the lore that if supersymmetry is not broken at tree-level, it remains 
unbroken in perturbation theory. The resolution is that the classical D-flat condition v u = Vd is 
perturbed, because of the isospin violation of the Yukawa couplings. The true flat direction is 
parameterized by: 

1 



16vr 2 



(4.4) 



3 J(m t 2 ) - 3 J(m 2 b ) - J(m 2 T ) 

One then finds that V c g is indeed up to terms of three-loop order; W ), contribute 
to V e g in the ratio 1 : —2 : 1. 



Scale dependence of the effective potential and running 

of A, v u and Vd 



An important consistency check on the effective potential formalism is that the value of Veg, 
being a physical quantity, should be independent of the arbitrary choice of renormalization scale 
Q. In the MSSM, the equation which expresses this is: 



Q -m v * 



du 



dv d 



cff 



0. 



(5.1) 



Here A represents all of the running DR' parameters of the theory (except the VEVs), and 

(S) (S) 

IhIiIhI are ^ ne anomalous dimensions of the Higgs scalar fields in Landau gauge. Note that 
these anomalous dimensions are gauge-dependent, and differ from the anomalous dimensions of 
the chiral superfields, because of gauge fixing. (See for example refs. [^8| and [p9[| .) Using the 



general formulas given in ref. [13], I obtain: 



(S) 



(S) 



1 



,(5,1) 



167T 2 
1 

16^2^. 



7^: + 

(.8,1) 



+ 



(167T 

1 



(16vr 



2)2 iHu ' 
(S,2) 



(5.2) 
(5.3) 
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where 



7 g,D = 3 Tr[Yt Yu ] - \g 2 - \g< 2 - (5.4) 
7 £ 2) = -9Tr[YtY u YtYj - 3Tr Y^ Y M Y^Y d ] + (16^ + ^5 , ' 2 )Tr[Y^ t Y u ] 

+3^ 4 + ^V 2 + f^ 4 ; (5.5) 



(5,i) 



3Tr[Y^Y d ] + TrYlYj - ^g 2 - \g' 2 - (5.6) 

7 g 2) = -9Tr[Y^Y,Y^Y d ] - 3Tr [Y+ Y u Yj Y d ] - 3Tr[Y|Y e Y|Y e ] 

9 ^9^ 
+ (16^ 2 - §^ ,2 )Tr[YtY d ] + 2^Tr[Y|Y e ] + 3g* + -^g' 2 + |</ 4 . (5.7) 



The DR VEVs run with renormahzation scale according to: 

= (5 ' 8) 
«^ = -4% (5.9) 

The DR' two-loop beta functions for all of the MSSM parameters can be found^ in ref. [[30 . 
The exception is the renormahzation group running of the DR' field-independent vacuum energy, 



which can be obtained from ref. 13 



/#> = 2(m 2 H J 2 + 2(m 2 H f + A\^m 2 Hu +m 2 H J + m 2 

+Tr[6mJm 2 3 + 2m 2 L m 2 L + 3m 2 m 2 + 3m 2 d m 2 d + m 2 e m 2 e ] 

-\ Ml \' - z\M2t - (s.ii) 

fif = 44^ 2 |M 1 | 4 + 36^|M 2 | 4 + | / i| 2 (24/|M 2 | 2 + 8^ 2 |M 1 | 2 ) 

+ (6/ + 2</ 2 ) {2\v\ 2 (m 2 Hu + m 2 Hi ) + «J 2 + (m 2 Hd ) 2 + 2|6| 2 + Tr[m 2 m 2 ]} 

2 16 
+(32<? 2 + 18/ + -( 7 ' 2 )Tr[m|m 2 3 ] + (16^ 2 + y ^ 2 )Tr[m 2 m 2 ] 

+ (16^ 2 + ^' 2 )Tr[m 2 m 2 ] + 4</ 2 Tr[m 2 m 2 ] 
-(12# 2 M 2 + \g' 2 M x )Vii - {l2g 2 M* + Ag' 2 M*)b/j,* 
-12Tr[Yt Y .J {(m 2 H J 2 + \^\ 2 (2m 2 Hu + m\ d ) + \b\ 2 } 

-(12Tr[YtYj+4Tr[Y e tY e ]) {(m 2 H f + \^\ 2 (2m 2 Hd + m 2 H J + \b\ 2 } 

^The parameters a u , a^, a e , and b in the present paper were represented by the symbols h u , h^, h e and B 
in ref. |30| . That reference also has an obvious overall minus sign error in the specification of the MSSM soft 
Lagrangian; £ should be —C in eqs. (4.2) and (4.3). The results given in that paper are actually in the DR 
scheme, although not explicitly identified as such at the time (see the "Note added" ) . 
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-12Tr[m 2 Q m 2 Q (Yt Y u + YjY d )] - 12Tr[m 2 m 2 Y u Yt] - 12Tr[m>M] 

-4Tr[m 2 m 2 Y^YJ - 4Tr[m 2 m 2 Y e Y+] 

-|/i| 2 {l2IV[ni 2 Q (Yl;Y u + YjY d )] + 12Tr[m 2 Y u Yt] 

+ 12Tr[m 2 Y d Yt] + 4Tr[m| Y|Y e ] + 4Tr[m 2 Y e Yt]} 

-12Tr[m 2 (ata, + aja^)] - 12Tr[m 2 a u at] - 12Tr[m 2 a d at] 

-4Tr[m|ata e ] -4Tr[m 2 a e at] - ^Tr^a^m 2 ^ + |/i| 2 ) 

-(121*^] +4TV[ata e ])(n4 + - 12Tr[Yt a J^* - 12Tr[4 Y„]//6 

-(12Tr[Y^a d ] + 4Tr[Yta e ])/ifc* - (12Tr[a+Y d ] + 4Tr[ a tY e ])/i*&. (5.12) 

Using these equations, I have checked that the effective potential found in section ||] is indeed 
renormalization scale-invariant at two-loop order. This demonstration is tedious and omitted. 

6 A numerical example 

I now present a quasi-realistic numerical example, to illustrate the results above. As a template 
model, I choose parameters at a renormalization scale Qq = 640 GeV: 

g > = 0.36, g = 0.65, g 3 = 1.06, y t = 0.90, y b = 0.13, y T = 0.10, 
M 1 = 150, M 2 = 280, M 3 = 800, a t = -600, a b = -150, a T = -40 GeV, 
m 2 Ql 2 = (780) 2 , ml i 2 = (740) 2 , m\ 2 = (735) 2 , m\ i 2 = (280) 2 , m 2 l2 = (200) 2 GeV 2 , 
(700) 2 , m 2 3 = (580) 2 , m\ = (725) 2 , mi 3 = (270) 2 , m% = (195) 2 GeV 2 , 
2 Hu = -(500) 2 , r< = (270) 2 GeV 2 (6.1) 



and 

(j, = 504.1811202 GeV; b = (184.2202586 GeV) 2 . (6.2) 

The last two values are engineered so that the minimum of the full 2-loop effective potential 
found in section § is: 

v u {Qo) = 172 GeV; v d {Q ) = 17.2 GeV. (6.3) 

One way to test the accuracy of the effective potential is by checking scale invariance. While 
I have done this analytically at two-loop order as described in section |j, in practice the neglected 
effects of higher order can be quite significant. To study this, I run the Lagrangian parameters of 
eqs. ( |6.1|) ,(|6.2| ) from the template scale Q to another scale Q, using the two- loop renormalization 
group equations of ref. fl30"| . The minimum of the effective potential at this new scale is found 



25 



176 
175 
174 
> 173 
172 
171 



17 



V"" 1 



2-loop V ff (Q) minimum 
2-loop running from Q =640 GeV I 



10.3 



10.2 - 



CO. 

U 101 



10.0 



00 200 300 400 500 600 700 800 900 1001 
Q [GeV] 



9: 



2-loop V rff (Q) minimum 

2-loop running from Q () =640 GeV 




?00 200 300 400 500 600 700 800 900 1000 
Q [GeV] 



Figure 2: Comparison of the values of v (Q) and tan (3(Q) as a test of renormalization scale invariance. 
The solid lines are values obtained directly from minimization of the two-loop effective potential at 
the scale Q, using Lagrangian parameters run from Qq = 640 GeV. The dashed lines are the result 



obtained by renormalization group running of the template VEVs of eq. (O) from Qq to Q using 
eqs. (|5.8[),(|5l)|). The discrepancy is due to effects at three loop order and beyond. 



numerically. The resulting VEVs are compared to the values obtained by running the values of 
eq. ( |6.3|) for v u and Vd obtained at Qo to Q using eqs. (|5\8| ),( ^9l) . The results of this comparison 
are shown in terms of the running quantities: 



v(Q) 



tan f3(Q) = v u /v d 



(6.4) 



in figure |2|. The good news is that the comparison of v(Q) obtained by the two methods shows 
agreement to better than 0.5% for a significant range of scales Q near the template scale Qq 
(which is approximately y/m^ mi 2 ) , and the comparison of tan/3 is good to better than 0.1%. 
However, near Q the slopes of v (Q) found by the two methods actually have the opposite 
sign! This surprising sensitivity to higher-loop effects is a consequence of the shallowness of the 
scalar potential along the direction v u /vd ~ tan/5. In contrast, higher-loop contributions have 
relatively much less effect on v u /vd, since that corresponds to a much steeper direction of the 
potential. 

An aside: the cusp-like feature found near Q = 463 GeV occurs because, at that scale, the 
tree-level squared mass of h° at the minimum of the two- loop effective potential goes through 0; 
it is positive for all larger Q. This leads to significant numerical effects because of the appearance 
of terms involving ln(m^ ) in the effective potential: 

1 



(mioYHmio/Q 2 ) + 



m 2 h0 



ci\n{m 2 h0 /Q 2 ) + c 2 ln 2 (m 2 h0 /Q 2 



+ .... 



(6.5) 
(6.6) 
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Figure 3: Comparisons of the values and scale dependence of the VEV parameters v and tan (3 found 
by minimizing the effective potential in different approximations. The label "min" refers to values 
obtained by direct minimization at Q, and "run" means obtained by running the template values of 
eq. (|6.3j ) from Qq. The solid lines follow from the full two- loop effective potential found in section 3 of 
the present paper. The dashed lines represent the one-loop approximation, and the dot-dashed lines 
the partial two-loop approximation of refs. 



where c\ and C2 have dimensions of (mass) 2 . Thus V e g is well-defined in the limit m\ — > 0. 
However, while the first derivatives of V^ 1 -* with respect to v u and Vd are finite, the first derivatives 
of have logarithmic and double-logarithmic divergences in that limit. 

A comparison of the full two-loop effective potential to previous approximations is shown 
in figure |3|. The graphs show the ratios v m i n /v run and tan/? min / tan/3 mn , where "min" means 
obtained by direct minimization at Q, and "run" means obtained by running the template 
values of v u , Vd in eq. ( |6.3| ) from Q to Q. This shows that the full two- loop effective potential 
indeed mitigates the scale dependence, and changes the value found for v compared to the 
previous state-of-the-art approximation in refs. by nearly 1.5% in this example. 

A different way to look at things, which may be closer to the situation we will face when 
confronted with experiment, is to view v u and Vd as input data rather than output parameters, 
and use the two minimization conditions of the effective potential to extract values for two 
other parameters. Since v u and Vd are especially sensitive to radiative corrections because of the 
shallowness of the potential, treating them as among the known quantities is advantageous. Here 
I will follow the commonplace procedure of treating b and fi as the unknowns to be solved for, 
with given values of v u , Vd and all the other Lagrangian parameters. (In reality, it seems clear that 
global fits to various observables will have to be conducted, since no more direct measurement 
of rn 2 Hu and m 2 Hd is possible.) In figure £|, I have graphed the ratios of /i m in//-trun and b m i n /b run , 
where "min" means that the effective potential at Q is required to be minimized, using values 
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Figure 4: Comparisons of the values and scale-dependence of the parameters // and b obtained by 
requiring the effective potential to be minimized. The label "min" refers to values obtained by running 
the Lagrangian parameters of eq. ( |6,lD and VEVs of eq. ( |6.3D to Q, and "run" means obtained by 
running the template values for fi, b of eq. ( |6. 2[ ) from Qq to Q. The solid lines follow from the full 
two-loop effective potential found in section 3. The dashed lines represent the one-loop approximation, 
and the dot-dashed lines the partial two- loop approximation of refs. 



obtained by running the Lagrangian parameters of eq. (|6.1| ) and VEVs of eq. (|6.3|) from Qq to Q, 
while "run" means those obtained by running the template values for /i, b of eq. ( |6.2|) from Qq 
to Q using two-loop renormalization group equations. In each case, the full two-loop results are 
compared to those obtained using one-loop and partial two-loop approximations for the effective 
potential. In this model, the scale-dependences of fi and b are less than a few hundredths of a 
percent, using the full two-loop potential over a wide range of scales Q. 



7 Outlook 

In this paper, I have presented the complete two-loop effective potential for the MSSM in the 
Landau gauge and the DR' scheme. 

The two-loop effective potential found here can, in principle, be renormalizat ion-group im- 
proved [J3TJ] to sum leading and sub-leading logarithms of ratios of different mass scales. However, 
the logarithmic contributions to the effective potential are typically not overwhelming compared 
to the non-logarithmic ones. Renormalization group improvement does give an improved scale 
dependence, because that is what it is designed to do. However, improved scale-independence 
does not always imply improved accuracy; it is a necessary but not sufficient criterion. There- 
fore, the efficacy of further renormalizat ion-group improvement of the MSSM effective potential 
on top of the full two-loop result is unclear to me. 
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If supersymmetry is part of the new physics associated with electroweak symmetry breaking, 
then these results will be part of a program to fit accurately experimental data to underlying 
Lagrangian parameters. Other parts of this program will require more precise calculations of 
the superpartner and Higgs scalar physical masses. Of particular importance is an improved 
calculation of the physical h° mass, which is well-known to be highly sensitive to radiative 
corrections. Unfortunately, m 2 h0 cannot simply be obtained by taking the second derivatives of 
the effective potential found here, because of important wave-function renormalization effects. 

The calculation described here can be extended to various non-minimal versions of supersym- 
metry, for example those with additional singlet fields. This can be done as a straightforward 
application of ref. flTJ| for a general theory. The fact that h° was not discovered at LEP may be 
taken as support for the importance of considering such models. 

This work was supported in part by the NSF grant number PHY-9970691. 
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